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Abstract 

Listing all the simple cycles (hereafter just called cycles) in a graph is a classical problem whose efficient 
solutions date back to the early 70s. For a graph with n vertices and m edges, containing r\ cycles, the 
\ best known solution in the literature is given by Johnson's algorithm [SIAM J. Computing, 1975] and takes 

0{{rj + l)(m + n)) time. This solution is surprisingly not optimal for undirected graphs: to the best of our 
^ ' knowledge, no theoretically faster solutions have been proposed in almost 40 years. 

We present the first optimal solution to list all the cycles in an undirected graph G, improving the time 
bound of Johnson's algorithm by a factor that can be 0(n 2 ). Specifically, let C(G) denote the set of all these 
cycles, and observe that C(G)| = to For a cycle c <E C(G), let c| denote the number of edges in c. Our 
algorithm requires 0(m + X)cec(G) l c l) time and is asymptotically optimal: indeed, fi(m) time is necessarily 
■ required to read G as input, and ^(X)ceC(G) l c l) time is necessarily required to list the output. 

We then describe an infinite family of dense graphs, in which each graph with n vertices and to edges 
contains r\ = J7(m) cycles c with |c| = O(l) edges. For each graph in this family, our algorithm requires 
0(ry + to) = 0(77) time, thus saving 0(m + n) — 0(n 2 ) time when compared to Johnson's algorithm. In 
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general for any graph, since \c\ < n, the cost of our algorithm never exceeds 0{m + (r/ + l)n) time. 

We also present the first optimal solution to list all the simple paths from s to f (shortly, si-paths) in 
an undirected graph G. Let V s t(G) denote the set of si-paths in G and, for an si-path it £ V s t(G), let |7r| 
be the number of edges in tt. Our algorithm lists all the si-paths in G optimally in 0(m + ^2 ne -p r G j |tt|) 
' time, observing that ^(X^gp t(G) time is necessarily required to list the output. 

While the basic approach is simple (see binary partition in point [3]), we use a number of non-trivial ideas 
to obtain our optimal algorithm for an undirected (connected) graph G as conceptually listed below. 

1. Prove the following reduction. If there exists an optimal algorithm to list the si-paths in G, there 
exists an optimal algorithm to list the cycles in G. This relates C(G) and V s t(G) for some choices s, t. 

2. Focus on listing the si-paths. Consider the decomposition of the graph into biconnected components 
(bccs), thus forming a tree T where two BCCs are adjacent in T iff they share an articulation point. 

■ Exploit (and prove) the property that if s and i belong to distinct BCCs, then (i) there is a unique 

sequence B s t of adjacent BCCs in T through which each si-path must necessarily pass, and (ii) each 
si-path is the concatenation of paths connecting the articulation points of these BCCs in B s t . 

3. Recursively list the si-paths in B Syt using the classical binary partition (i.e. given an edge e in G, list 
all the cycles containing e, and then all the cycles not containing e): now it suffices to work on the first 
BCC in B s j, and efficiently maintain it when deleting an edge e, as required by the binary partition. 

4. Use a notion of certificate to avoid recursive calls (in the binary partition) that do not list new si-paths. 
This certificate is maintained dynamically as a data structure representing the first BCC in B s t ^ which 
guarantees that there exists at least one new solution in the current B s t . 

5. Consider the binary recursion tree corresponding to the binary partition. Divide this tree into spines: 
a spine corresponds to the recursive calls generated by the edges e belonging to the same adjacency 
list in B s .t- The amortized cost for each listed si-path 7r is 0(|7r|) if there is a guarantee that the 
amortized cost in each spine S is O(fi), where /1 is a lower bound on the number of si-paths that 
will be listed from the recursive calls belonging to S. The (unknown) parameter (i different for each 
spine S, and the corresponding cost 0(/i), will drive the design of the proposed algorithms. 



1 Introduction 



Listing all the simple cycles (hereafter just called cycles) in a graph is a classical problem whose 
efficient solutions date back to the early 70s. For a graph with n vertices and m edges, containing 
rj cycles, the best known solution in the literature is given by Johnson's algorithm [7J and takes 
0((n + l)(m + n)) time. This solution is surprisingly not optimal for undirected graphs: to the 
best of our knowledge, no theoretically faster solutions have been proposed in almost 40 years. 

Our results. We present the first optimal solution to list all the cycles in an undirected graph G, 
improving the time bound of Johnson's algorithm by a factor that can be 0(n 2 ). Specifically, let 
C(G) denote the set of all these cycles, and observe that |C(G)| = rj. For a cycle c € C(G), 
let |c| denote the number of edges in c. Our algorithm requires 0{m + X^ceC(G) l c l) time and 
is asymptotically optimal: indeed, f2(m) time is necessarily required to read G as input, and 
^(SceC(G) l c D time is necessarily required to list the output. Since |c| < n, the cost of our algorithm 
never exceeds 0(m + (77 + l)n) time. 

Along the same lines, we also present the first optimal solution to list all the simple paths from 
s to t (shortly, si-paths) in an undirected graph G. Let V s t{G) denote the set of si-paths in G and, 
for an si-path tt £ V s t(G), let \tt\ be the number of edges in it. Our algorithm lists all the si-paths 
in G optimally in 0{m + ^neVstiG) M) time, observing that ^(X^e-p^G) I 71 "!) time is necessarily 
required to list the output. We prove the following reduction to relate C{G) and V s t{G) for some 
suitable choices of vertices s,i: If there exists an optimal algorithm to list the si-paths in G, then 
there exists an optimal algorithm to list the cycles in G. Hence, we can focus on listing si-paths. 

History of the problem. The classical problem of listing all the cycles of a graph has been 
extensively studied for its many applications in several fields, ranging from the mechanical analysis 
of chemical structures [16] to the design and analysis of reliable communication networks, and the 
graph isomorphism problem |21j . There is a vast body of work, and the majority of the algorithms 
listing all the cycles can be divided into the following three classes (see [HE] for excellent surveys). 

Search space algorithms. According to this approach, cycles are looked for in an appropriate 
search space. In the case of undirected graphs, the cycle vector space [3] turned out to be the most 
promising choice: from a basis for this space, all vectors are computed and it is tested whether 
they are a cycle. Since the algorithm introduced in [21] . many algorithms have been proposed: 
however, the complexity of these algorithms turns out to be exponential in the dimension of the 
vector space, and thus in n. For the special case of planar graphs, in [17] the author was able to 
design an algorithm listing all the cycles in 0((rj + l)n) time. 

Backtrack algorithms. According to this approach, all paths are generated by backtrack and, for 
each path, it is tested whether it is a cycle. One of the first algorithms based on this approach is the 
one proposed in [20], which is however exponential in rj. By adding a simple pruning strategy, this 
algorithm has been successively modified in [19]: it lists all the cycles in 0(nm(r]+l)) time. Further 
improvements were proposed in [7J, [18], and [13], leading to 0((rj + l)(m + n))-time algorithms 
that work for both directed and undirected graphs. 

Algorithms using the powers of the adjacency matrix. This approach uses the so-called variable 
adjacency matrix, that is, the formal sum of edges joining two vertices. A non-zero element of 
the p-th. power of this matrix is the sum of all walks of length p: hence, to compute all cycles, 
we compute the nth power of the variable adjacency matrix. This approach is not very efficient 
because of the non-simple walks. All algorithms based on this approach (e.g. [12] and [23]) basically 
differ only on the way they avoid to consider walks that are neither paths nor cycles. 

Almost 40 years after Johnson's algorithm [7J, the problem of efficiently listing all cycles of a 
graph is still an active area of research (e.g. [21 El [TO"! [T4"l [22| [15]). New application areas have 
emerged in the last decade, such as bioinformatics: for example, two algorithms for this problem 
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have been proposed in [8] and [9] while studying biological interaction graphs. Nevertheless, no 
significant improvement has been obtained from the theory standpoint: in particular, Johnson's 
algorithm is still the theoretically most efficient. 

Hard graphs for Johnson's algorithm. We now describe an infinite family of dense undi- 
rected graphs, in which each graph with n vertices and m edges contains n = Q(rn) cycles c with 
|c| = O(l) edges. Suppose w.l.o.g. that n is a multiple of 3, and consider a tripartite complete 
graph = (V% U V% U Vs,E), where V\, V2, V3 are pairwise disjoint sets of vertices of size n/3 

each, and E is the set of edges thus formed: for each choice of x\ G V\,xi G ^,£3 G V3, edges 
(x\,X2), (x2,xs), and (xi,xs) belong to E; and no edges exist that connect any two vertices within 
Vi, for i = 1,2,3. Note that each choice of x\ G V±,X2 G ^,£3 G V3 in K^, 3 gives rise to a 
distinct cycle c = (x\,X2), (22, £3), (xs, xi) of length |c| = 3, and there are no other (simple) cycles. 
Thus, there are rj = (n/3) 3 cycles in the graph -f^ 3 / 3 , where m = 3(n/3) 2 . For each graph K^ 3 in 
this family, our algorithm requires 0(n + m) = ©(77) = 0(n 3 ) time to list the cycles, thus saving 
0(m + n) = 0(n 2 ) time when compared to Johnson's algorithm that takes 0(n 5 ) in this case. 

Note that the analysis of the time complexity of Johnson's algorithm is not pessimistic and can- 
not match the one of our algorithm for listing cycles. For example, consider the sparse "diamond" 
graph D n = (V, E) in Fig.[T]with n = 2k+3 vertices in V = {a, b, c, v\, . . . , Vk, u±, . . . , u^}. There are 
m = Q(n) edges in E = {(a, c), (a, Vi), (vi, b), (b, U{), (ui, c), for 1 < i < k}, and three kinds of (sim- 
ple) cycles: (1) (a,Vi), (vi,b), (b,Uj), (uj,c), (c,a) for 1 < i, j < k; (2) (a,Vi), (vi,b), (b,Vj), (vj,a) for 

1 < i < j < k; (3) (b,Ui), (ui,c), (c,Uj), (uj,b) for 1 < i < j < k, totalizing rj = 0(n 2 ) cycles. Our 
algorithm takes 0(n + k 2 ) = @(rj) = 0(n 2 ) time to list these cycles. On the other hand, Johnson's 
algorithm takes 0(n 3 ) time, and the discovery of the 0(n 2 ) cycles in (1) costs 0(/c) = 0(n) time 
each: the backtracking procedure in Johnson's algorithm starting at a, and passing through Vi, b 
and Uj for some i, j, arrives at c: at that point, it explores all the vertices ui (I 7^ i) even if they do 
not lead to cycles when coupled with a, Vi, b, Uj, and c. 

2 Overview and Main Ideas 
Preliminaries 

Let G = (V, E) be an undirected connected graph with n = \V\ vertices and m = \E\ edges, without 
self-loops or parallel edges. For a vertex u G V, we denote by N(u) the neighborhood of u and by 
d(u) = \N(u) \ its degree. G[V] denotes the subgraph induced by V C V, and G — u is the induced 
subgraph G[V \{u}] for u G V. Likewise for edge e G E, we adopt the notation G — e = (V, E\{e}). 

Paths are simple in G by definition: we refer to a path tt by its natural sequence of vertices 
or edges. A path ir from s to t, or st-path, is denoted by ir = s ~^ t. Additionally, V{G) is the 
set of all paths in G and V s .t(G) is the set of all si-paths in G. When s = t we have cycles, and 
C(G) denotes the set of all cycles in G. We denote the number of edges in a path tt by \ir\ and in 
a cycle c by |c|. In this paper, we consider the following problems. 

Problem 1 (Listing st-Paths) Given a graph G = (V, E) and two distinct vertices s,t G V , 
output all the paths tt G V s ,t{G). 

Problem 2 (Listing Cycles) Given a graph G = (V,E), output all the cycles c G C{G). 

Our algorithms assume without loss of generality that the input graph G is connected, hence 
m > re—1, and use the decomposition of G into biconnected components. Recall that an articulation 
point (or cut- vertex) is a vertex u £ V such that the number of connected components in G increases 
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Figure 1: Diamond graph. 



Figure 2: Block tree of G with bead string B s ^ in gray. 



when u is removed. G is biconnected if it has no articulation points. Otherwise, G can always be 
decomposed into a tree of biconnected components, called the block tree, where each biconnected 
component is a maximal biconnected subgraph of G (see Fig. [2]), and two biconnected components 
are adjacent if and only if they share an articulation point. 

2.1 Reduction to si-paths 

We now show that listing cycles reduces to listing si-paths while preserving the optimal complexity. 

Lemma 3 Given an algorithm that solves ProblemUlin optimal 0(m + Yl^eVs t (G) M) time, there 
exists an algorithm that solves Problem\^in optimal 0(m + ^ ce c(G) l c l) time. 

Proof: Compute the biconnected components of G and keep them in a list L. Each (simple) cycle 
is contained in one of the biconnected components and therefore we can treat each biconnected 
component individually as follows. While L is not empty, extract a biconnected component B = 
(Vb,Eb) from L and repeat the following three steps: (i) compute a DFS traversal of B and take 
any back edge b = (s, t) in B; (ii) list all si-paths in B — b, i.e. the cycles in B that include edge b; 
(Hi) remove edge b from B, compute the new biconnected components thus created by removing 
edge b, and append them to L. When L becomes empty, all the cycles in G have been listed. 

Creating L takes 0(m) time. For every B E L, steps (i) and (Hi) take 0(\Eb\) time. Note that 
step (ii) always outputs distinct cycles in B (i.e. si-paths in B—b) in 0(\EB\+Y^ n& -p a t(B-b) time. 
However, B — b is then decomposed into biconnected components whose edges are traversed again. 
We can pay for the latter cost: for any edge e ^ b in a biconnected component B, there is always a 
cycle in B that contains both b and e (i.e. it is an si-path in B — b), hence Ylir£V a t (B-b) M dominates 
the term \E B \, i.e. ^2 weVs t ( B -b) M = ^(\ e b\)- Therefore steps take 0(Y,n£Vs,t(B-b) H) 

time. When L becomes empty, the whole task has taken 0(m + X^ceC(G) l c l) time. □ 

2.2 Decomposition into biconnected components 

We now focus on listing si-paths (Problem [1]) . We use the decomposition of G into a block tree of 
biconnected components. Given vertices s,i, define its bead string, denoted by B s j, as the unique 
sequence of one or more adjacent biconnected components (the beads) in the block tree, such that 
the first one contains s and the last one contains i (see Fig. [2]): these biconnected components are 
connected through articulation points, which must belong to all the paths to be listed. 

Lemma 4 All the st-paths in V Si t(G) are contained in the induced subgraph G[B s j] for the bead 
string B s j. Moreover, all the articulation points in G[B S ^} are traversed by each of these paths. 

Proof: Consider an edge e = (u,v) in G such that u E B s ^ and v ^ B S)t - Since the biconnected 
components of a graph form a tree and the bead string B s j is a path in this tree, there are no 
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paths v w in G — e for any w G B s j because the biconnected components in G are maximal 
and there would be a larger one (a contradiction). Moreover, let B%, B2, ■ ■ ■ , B r be the biconnected 
components composing B s j, where s € B\ and t G B r . If there is only one biconnected component 
in the path (i.e. r = 1), there are no articulation points in B s j- Otherwise, all of the r — 1 
articulation points in B s t are traversed by each path n € V s ,t(G): indeed, the articulation point 
between adjacent biconnected components Bi and -Bj+i is their only vertex in common and there 
are no edges linking Bi and flj+i. □ 

We thus restrict the problem of listing the paths in V s ,t{G) to the induced subgraph G[B St t], 
conceptually isolating it from the rest of G. For the sake of description, we will use interchangeably 
B s j and G[B s j] in the rest of the paper. 

2.3 Binary partition scheme 

We list the set of si-paths in B s j, denoted by V St t(B s j), by applying the binary partition method 
(where V s ,t(G) = V s ,t(Bs,t) by LemmaH]): We choose an edge e = (s,v) incident to s and then list 
all the si-paths that include e and then all the si-paths that do not include e. Since we delete some 
vertices and some edges during the recursive calls, we proceed as follows. 

Invariant: At a generic recursive step on vertex u (initially, u := s), let ir s = s ~» u be the 
path discovered so far (initially, ir s is empty {}). Let B u j be the current bead string (initially, 
B U)t := B St t). More precisely, B U)t is defined as follows: (i) remove from B s t all the nodes in ir s but 
u, and the edges incident to u and discarded so far; (ii) recompute the block tree on the resulting 
graph; (Hi) B U)t is the unique bead string that connects u to t in the recomputed block tree. 

Base case: When u = t, output the st-path ir s computed so far. 

Recursive rule: Let V(tt s ,u, B u j) denote the set of si-paths to be listed by the current recursive 
call. Then, it is the union of the following two disjoint sets, for an edge e = (u, v) incident to u: 

• [left branching] the st-paths in "P(7r s • e,v,B v j) that use e, where B v j is the unique bead 
string connecting v to t in the block tree resulting from the deletion of vertex u from B u t] 

• [right branching] the st-paths in V(ir s , u, B' u t ) that do not use e, where B' ut is the unique 
bead string connecting u to t in the block tree resulting from the deletion of edge e from B u ^- 

Hence, V s ,t(B s ,t) (and so V s ,t{G)) can be computed by invoking V({}, s, B Sj t)- The correctness and 
completeness of the above approach is discussed in Section 12.41 

At this point, it should be clear why we introduce the notion of bead strings in the binary 
partition. The existence of the partial path ir s and the bead string B u j guarantees that there 
surely exists at least one si-path. But there are two sides of the coin when using B u t- 

(1) One advantage is that we can avoid useless recursive calls: If vertex u has only one incident 
edge e, we just perform the left branching; otherwise, we can safely perform both the left and right 
branching since the first bead in B u j, is always a biconnected component by definition (and so there 
exist both an si-path that traverses e and one that does not traverse e). 

(2) The other side of the coin is that we have to maintain the bead string B u j as B v $ in the 
left branching and as B' u t in the right branching by Lemma HI Note that these bead strings are 
surely non-empty since B u j is non-empty by induction (we only perform either left or left /right 
branching when there are solutions by point (!))■ 

To efficiently address point (2), we need to introduce the notion of certificate as described next. 
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Figure 3: Example certificate C 



Algorithm 1 list_paths s t (ir s , u, C) 

1: if u = t then 
2: output(-7r s ) 
3: return 

4: end if 

5: e = (u,v) := choose(C, u) 

6: if e is back edge then 
7: I := right_update(C, e) 
8: list_paths s t (7r s , u, C) 
9: restore(C, /) 

10: end if 

11: I := left_update(C, e) 

12: list_paths s t (7r s • («, i>), t>, C) 

13: restore(C, I) 



2.4 Introducing the certificate 

Given the bead string B Ujt , we call the head of B Ut t, denoted by H u , the first biconnected component 
in B Ut t, where u € H u . Consider a DFS tree of B U)t rooted at u that changes along with B u j, and 
classify the edges in B u> t as tree edges or back edges (no cross edges since the graph is undirected). 

To maintain B Ut t (and so H u ) during the recursive calls, we introduce a certificate C (see Fig. [3] 
for an example): It is a suitable data structure that uses the above classification of the edges in 
B u j, and supports the following operations, required by the binary partition scheme. 

• choose(C, u): returns an edge e = (u, v) with v € H u such that tt s ■ (u, v) ■ u ~~> t is an st-path 
such that u ~^ t is inside B u ^. Note that e always exists since H u is biconnected. Also, the 
chosen v is the last one in DFS order among the neighbors of it: in this way, the (only) tree 
edge e is returned when there are no back edges leaving from uQ 

• lef t_update(C, e): for the given e = (u,v), it obtains B v> t from B U) t as discussed in Sec- 
tion [231 This implies updating also H u , C, and the block tree, since the recursion continues 
on v. It returns bookkeeping information / for what is updated, so that it is possible to revert 
to B u j, H u , C, and the block tree, to their status before this operation. 

• right_update(C, e): for the given e = (u,v), it obtains B' ut from B u t as discussed in Sec- 
tion 12.31 which implies updating also H u , C, and the block tree. It returns bookkeeping 
information / as in the case of lef t_update(C, e). 

• restore(C, I): reverts the bead string to B u t, the head H u , the certificate C, and the block 
tree, to their status before operation I := lef t_update(C, e) or / := right_update(C, e) was 
issued (in the same recursive call). 

Note that a notion of certificate in listing problems has been introduced in [4J, but it cannot 
be directly applied to our case due to the different nature of the problems and our use of more 
complex structures such as biconnected components. Using our certificate and its operations, we 
can now formalize the binary partition and its recursive calls V(ir s ,u, B u j) described in Section [231 
as Algorithm [H where B u t is replaced by its certificate C. 

As it will be clear in Sections [3] and [4] this order facilitates the analysis and the implementation of the certificate. 
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The base case (u = t) corresponds to lines 1-4 of Algorithm [TJ During recursion, the left 
branching corresponds to lines 5 and 11-13, while the right branching to lines 5-10. Note that we 
perform only the left branching when there is only one incident edge in u, which is a tree edge by 
definition of choose. Also, lines 9 and 13 are needed to restore the parameters to their values when 
returning from the recursive calls. The proof of the following lemma is in the Appendix. 

Lemma 5 Algorithm^ correctly lists all the st-paths in V s ,t{G). 

A natural question is what is the complexity: we should account for the cost of maintaining C 
and for the cost of the recursive calls of Algorithm [Tj Since we cannot always maintain the certificate 
in 0(1) time, the ideal situation for attaining an optimal cost is taking 0(/u) time if at least \x st- 
paths are listed in the current call (and its nested calls). An obstacle to this ideal situation is that 
we cannot estimate \i efficiently and cannot design Algorithm Q] so that it takes O(p) adaptively. 
We circumvent this by using a different cost scheme in Section 12.51 that is based on the recursion 
tree induced by AlgorithmQ] Section|3]and Appendix|B]are devoted to the efficient implementation 
of the above certificate operations according to the cost scheme that we discuss next. 

2.5 Recursion tree and cost amortization 

We now define how to distribute the costs among the several recursive calls of Algorithm [TJ so that 
optimality is achieved. Consider a generic execution on the bead string B Ui t- We can trace this 
execution by using a binary recursion tree R. The nodes of R are labeled by the arguments on 
which Algorithm Q] is run: specifically, we denote a node in R by the triple x = (tt s ,u,C) iff it 
represents the call with arguments ir s , u, and cH The left branching is represented by the left 
child, and the right branching (if any) by the right child of the current node. 

Lemma 6 The binary recursion tree R for B u j has the following properties: 

1. There is a one-to-one correspondence between the paths in V s ,t(B Ut t) and the leaves in the 
recursion tree rooted at node (ir s ,u,C). 

2. Consider any leaf and its corresponding st-path ir: there are \tt\ left branches in the corre- 
sponding root-to-leaf trace. 

3. Consider the instruction e := choose(C, u) in Algorithm^ unary (i.e. single-child) nodes 
correspond to left branches (e is a tree edge) while binary nodes correspond to left and right 
branches (e is a back edge). 

4- The number of binary nodes is [P s j(B u t)\ — 1. 

We define a spine of R to be a subset of -R's nodes linked as follows: the first node is a node x 
that is either the left child of its parent or the root of R, and the other nodes are those reachable 
from x by right branching in R. Let x = (ir s ,u,C) be the first node in a spine S. The nodes in S 
correspond to the edges that are incident to vertex u in B u ^: hence their number equals the degree 
d(u) of u in B U) t, and the deepest (last) node in S is always a tree edge in B u j while the others 
are back edges (in reverse DFS order). Summing up, R can be seen as composed by spines, unary 
nodes, and leaves, where each spine has a unary node as deepest node. This gives a global pictures 
of R that we now exploit for the analysis. 

2 For clarity, we use "nodes" when referring to R and "vertices" when referring to B Utt - 
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We define the compact head, denoted by Hx = (Vx,Ex), as the (multi)graph obtained by 
compacting the maximal chains of degree-2 vertices, except u, t, and the vertices that are the 
leaves of its DFS tree rooted at u. 

The rationale behind the above definition is that the costs defined in terms of Hx amortize 
well, as the size of Hx and the number of si-paths in the subtree of R rooted at node x = (ir s ,u, C) 
are intimately related (see Lemma [10] in Section [3]) while this is not necessarily true for H u . 

We now define the following abstract cost for spines, unary nodes, and leaves of R, for a 
sufficiently large constant cq > 0, that Algorithm Q] must fulfill: 

{Co if v is a unary node 

co|vr| if r is a leaf corresponding to path tt (1) 

Co (| Vx | + \Ex\) if t is a spine with compact head Hx 

Lemma 7 The sum of the costs in the nodes of the recursion tree X^re_R^( r ) = ^C^nreVs t{B u t ) I 71 "!)' 

Section [3] contains the proof of Lemma[7]and related properties. Setting u := s, we obtain that 
the cost in Lemma [7] is optimal, by Lemma [H 

Theorem 8 Algorithm^ solves problem Problem^ in optimal 0{m + ^2 7T£ -p s t ng\ M) time. 

By Lemma El we obtain an optimal result for listing cycles. 
Theorem 9 Problem^ can be optimally solved in 0(m + X^ceC(G) l c l) time. 

3 Amortization strategy 

We devote this section to prove LemmaEl Let us split the sum in Eq. (pQ) in three parts, and bound 
each part individually, as 

£>(r)< £ T(r)+Y,T(r)+ £ T(r). (2) 

rgiJ r: unary r: leaf r: spine 

We have that ^r- unary -^( r ) = ^(Z^TrGPs t {G) I 71 "!)' smce there are \V s ,t(G)\ leaves, and the root- 
to-leaf trace leading to the leaf for ir contains at most \tt\ unary nodes by Lemma [6j where each 
unary node has cost 0(1) by Eq. ([T]). 

Also, Zlr:ieaf^( r ) = ^(DttsT's t (G) I 71 "!)' snice tne l e& f r f° r 7T has cost O ( 1 7T | ) by Eq. (P). 

It remains to bound ^ r spine T(r). By Eq. ([T]), we can rewrite this cost as Y2h x cod^x + \Ex\), 
where the sum ranges over the compacted heads Hx associated with the spines r. We use the 
following lemma to provide a lower bound on the number of st-paths descending from r. 

Lemma 10 Given a spine r, and its bead string B u t with head H u , there are at least \Ex \ — \Vx \ + 1 
st-paths in G that have prefix ir s = s ~» u and suffix u ~» t internal to B u t, where the compacted 
head is Hx = (Vx,Ex). 

Proof: Hx is biconnected. In any biconnected graph B = (Vb,Eb) there are at least \Eb \ — \ Vb\ + 1 
xy-paths for any x,y € Vb- Find a ear decomposition [3] of B and consider the process of forming B 
by adding ears one at the time, starting from a single cycle including x and y. Initially \Vb \ = \Eb\ 
and there are 2 xy-paths. Each new ear forms a path connecting two vertices that are part of a 
xy-path, increasing the number of paths by at least 1. If the ear has k edges, its addition increases 
V by k — 1, E by k, and the number of xy-paths by at least 1. The result follows by induction. □ 
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The implication of Lemma PTOl is that there are at least \Ex\ — \Vx\ + 1 leaves descending from 
the given spine r. Hence, we can charge to each of them a cost of i^-pi^^j ■ Lemma HU allows 
us to prove that the latter cost is 0(1) when H u is different from a single edge or a cycle. (If H u 
is a single edge or a cycle, Hx is a single or double edge, and the cost is trivially a constant.) 

Lemma 11 For a compacted head Hx = (Vx,Ex), its density is > ^j. 

Specifically, let a = ^ and write a = 1 + 2/(3 for a constant /?: we have that \Ex \ + \Vx\ = 
(\E X \ ~ \V X \) + 2\V X \ < (\E X \ - \Vx\)+P(\E x \ - \V X \) = *%(\E X \ ~ \V X \). Thus, we can charge 
each leaf with a cost of Tg^pj^^^ < c oi~r = ^ (-*-)■ This motivates the definition of Hx, since 
Lemma [TT1 does not necessarily hold for the head H u (due to the unary nodes in its DFS tree). 

One last step to bound Y1h x c o(l^d + |-^x|) : as noted before, a root-to-leaf trace for the string 
storing ir has |7r| left branches by Lemma [6j and as many spines, each spine charging co^j = 0(1) 
to the leaf at hand. This means that each of the \V s ,t(G) \ leaves is charged for a cost of 0(|-7r|), thus 
bounding the sum as Er Sp ine T ( r ) = Eff x c (jVx| + \E X \) = Od^-gp^c) M). This completes 
the proof of Lemma [71 As a corollary, we obtain the following result. 

Lemma 12 The recursion tree R with cost as in Eq. ([I]) induces an 0(\tt\) amortized cost for each 
st-path ir. 

4 Certificate implementation and maintenance 

The certificate C associated with a node {tt s ,u,C) in the recursion tree is a compacted and aug- 
mented DFS tree of bead string B u t, rooted at vertex u. The DFS tree changes over time along 
with B u> t, and is maintained in such a way that t is in the leftmost path of the tree. We compact 
the DFS tree by contracting the vertices that have degree 2, except u, t, and the leaves (the latter 
surely have incident back edges). Maintaining this compacted representation is not a difficult data- 
structure problem. From now on we can assume w.l.o.g. that C is an augmented DFS tree rooted 
at u where internal nodes of the DFS tree have degree > 3, and each vertex v has associated: 

1. A doubly-linked list lb(v) of back edges linking v to its descendants w sorted by DFS order. 

2. A doubly-linked list ab(v) of back edges linking v to its ancestors w sorted by DFS order. 

3. An integer j(v), such that if v is an ancestor of w then j(v) < 7(10). 

4. The smallest j(w) over all neighbors w of v in C, excluding the parent, denoted by lowpoint(v). 

Given two vertices v,w € C such that v is the parent of w, we can efficiently test if any of the 
children of w is in the same biconnected component of v, i.e. lowpoint(w) < j(v). (Note that we 
adopt a variant of lowpoint using 7(1;) in place of depth(v): it has the same effect whereas using 
"f(v) is preferable since it is easier to dynamically maintain.) 

Lemma 13 The certificate associated with the root of the recursion can be computed in 0(m) time. 

Proof: In order to set t to be in the leftmost path, we perform a DFS traversal of graph G starting 
from s and stop when we reach vertex t. We then compute the DFS tree, traversing the path 
s ~» t first. When visiting vertex v, we set j(v) to depth of v in the DFS. Before going up on the 
traversal, we compute the lowpoints using the lowpoints of the children. Let z be the parent of v. 
If lowpoint (v) > lowpoint (z) and w is not in the leftmost path in the DFS, we cut the subtree of 
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v as it does not belong to B Syt . When finding a back edge e = (v, w), if w is a descendant of v we 
append e to both lb{v) and ab(w); else we append e to both afe(f) and lb(w). This maintains the 
DFS order in the back edge lists. This procedure takes at most two DFS traversals in 0(m) time. 
This DFS tree can be compacted in the same time bound. □ 

Lemma 14 Operation choose(C, u) can be implemented in O(l) time. 

Proof: If the list lb(v) is empty, return the tree edge e = (u,v) linking u to its only child v (there 
are no other children). Else, return the last edge in lb{v). □ 

We analyze the cost of updating and restoring the certificate C. We can reuse parts of C, 
namely, those corresponding to the vertices that are not in the compacted head Hx = (Yx,Ex) 
as defined in Section 1231 We prove that, given a unary node u and its tree edge e = (u,v), the 
subtree of v in C can be easily made a certificate for the left branch of the recursion. 

Lemma 15 On a unary node, lef t_update(C, e) takes 0(1) time. 

Proof: Take edge e = (u,v). Remove edge e and set v as the root of the certificate. Since e is the 
only edge incident in v, the subtree v is still a DFS tree. Cut the list of children of v keeping only 
the first child. (The other children are no longer in the bead string and become part of J.) There 
is no need to update 7(f)- D 

We now devote the rest of this section to show how to efficiently maintain C on a spine. Consider 
removing a back edge e from u: the compacted head Hx = (Vx,Ex) of the bead string can be 
divided into smaller biconnected components. Many of those can be excluded from the certificate 
(i.e. they are no longer in the new bead string, and so they are bookkept in I) and additionally 
we have to update the lowpoints that change. We prove that this operation can be performed in 
0(|Vx|) total time on a spine of the recursion tree. 

Lemma 16 The total cost of all the operations right_update(C, e) in a spine is 0(|Vx|) time. 

Proof: In the right branches along a spine, we remove all back edges in lb{u). This is done by 
starting from the last edge in lb(u), i.e. proceeding in reverse DFS order. For back edge b% = (zi,u), 
we traverse the vertices in the path from Z{ towards the root u, as these are the only lowpoints that 
can change. While moving upwards on the tree, on each vertex w, we update lowpoint (w). This is 
done by taking the endpoint y of the first edge in ab(w) (the back edge that goes the topmost in 
the tree) and choosing the minimum between j(y) and the lowpoint of each child of w. We stop 
when the updated lowpoint (w) = j(u) since it implies that the lowpoint of the vertex can not be 
further reduced. Note that we stop before u, except when removing the last back edge in lb(u). 

To prune the branches of the DFS tree that are no longer in B u j, consider again each vertex w 
in the path from Z{ towards the root u and its parent y. It is possible to test if w is an articulation 
point by checking if the updated lowpoint(w) > j(y). If that is the case and w is not in the leftmost 
path of the DFS, it implies that w ^ B Ut t, and therefore we cut the subtree of w and keep it in I 
to restore later. We use the same halting criterion as in the previous paragraph. 

The cost of removing all back edges in the spine is 0(|Vx|): there are 0(|Vx|) tree edges and, 
in the paths from Z{ to u, we do not traverse the same tree edge twice since the process described 
stops at the first common ancestor of endpoints of back edges \. Additionally, we take 0(1) time 
to cut a subtree of an articulation point in the DFS tree. □ 
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To compute lef t_update(C, e) in the binary nodes of a spine, we use the fact that in every left 
branching from that spine, the graph is the same (in a spine we only remove edges incident to u 
and on a left branch from the spine we remove the node u) and therefore its block tree is also the 
same. However, the certificates on these nodes are not the same, as they are rooted at different 
vertices. By using the DFS order of the edges, we are able to traverse each edge in Hx only a 
constant number of times in the spine. 

Lemma 17 The total cost of all operations lef t_update(C, e) in a spine is amortized 0(\Ex\)- 

Proof: Let t' be the last vertex in the path u ~^ t s.t. t' E Vx- Since t' is an articulation point, 
the subtree of the DFS tree rooted in t' is maintained in the case of removal of vertex u. Therefore 
the only modifications of the DFS tree occur in the compacted head Hx of B u ^. Let us compute 
the certificate Cf this is the certificate of the left branch of the ith node of the spine where we 
augment the path with the back edge bi = (Zi,u) of lb(u) in the order defined by choose(C, u). 

For the case of Ci, we remove u and rebuild the certificate starting form %\ (the last edge in 
lb(u)) using the algorithm from Lemma [131 restricted to Hx and using t' as target and 7(t') as a 
baseline to 7 (instead of the depth). This takes 0(\Ex\)- 

For the general case of Cj with i > 1 we also rebuild (part) of the certificate starting from z% 
using the procedure from Lemma [13] but we use information gathered in Cj_i to avoid exploring 
useless branches of the DFS tree. The key point is that, when we reach the first bead in common 
to both B Zu t and B Zil j, we only explore edges internal to this bead. If an edge e leaving the bead 
leads to t, we can reuse a subtree of Cj_i. If e does not lead to t, then it has already been explored 
(and cut) in Q_i and there is no need to explore it again since it will be discarded. Given the order 
we take bi, each bead is not added more than once, and the total cost over the spine is 0{\Ex\)- 

Nevertheless, the internal edges E' x of the first bead in common between B Zi j and B Zi _ x ^ can 
be explored several times during this procedurell We can charge the cost 0(\E' X \) of exploring 
those edges to another node in the recursion tree, since this common bead is the head of at least 
one certificate in the recursion subtree of the left child of the ith node of the spine. Specifically, 
we charge the first node in the leftmost path of the ith node of the spine that has exactly the 
edges E' x as head of its bead string: (i) if \E' X \ < 1 it corresponds to a unary node or a leaf in the 
recursion tree and therefore we can charge it with O(l) cost; (ii) otherwise it corresponds to a first 
node of a spine and therefore we can also charge it with 0(\E' X \). We use this charging scheme 
when i 7^ 1 and the cost is always charged in the leftmost recursion path of ith node of the spine. 
Consequently, we never charge a node in the recursion tree more than once. □ 

Lemma 18 On each node of the recursion tree, restore(C, /) takes time proportional to the size 
of the modifications kept in I. 

From Lemmas 1141 and 1161118} it follows that on a spine of the recursion tree we have the costs: 
choose(u) on each node which is bounded by 0(|Vx|) as there are at most \Vx\ back edges in u; 
right_update(C, e), restore(C, I) take 0(|Vx|) time; lef t_update(C, e) and restore(C, /) are 
charged 0{\ Vx | + |-Sx|) time. We thus have the following result, completing the proof of Theorem[5J 

Lemma 19 Algorithm^ can be implemented with a cost fulfilling Eq. ([1]), thus it takes total 0(m + 
Erei?. T ( r )) = °( m + Y.^VsABu,t) M) Ume - 



3 Consider the case where Zi,...,Zj are all in the same bead after the removal of u. The bead strings are the same, 
but the roots Zi, . . . , Zj are different, so we have to compute the corresponding DFS of the first component \j — i\ 
times. 
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A Omitted Proofs 



A.l Lemma [5] 

Proof: For a given vertex u the function choose(C, u) returns an edge e incident to u. We maintain 
the invariant that ir s is a path s u, since at the point of the recursive call in line [121 (i) is 
connected as we append edge (it, u) to tt s and; (ii) it is simple as vertex u is removed from the 
graph G in the call to lef t_update(C, e) in line llll In the case of recursive call in line [8] the 
invariant is trivially maintained as n s does not change. The algorithm only outputs si-paths since 
7r s is a s u path and u = t when the algorithm outputs, in line [2 

All the paths with prefix tt s that do not use e are listed by the recursive call in line EJ This is 
done by removing e from the graph in line [7] and thus no path can include e. The paths that use e 
are listed in line [12] since in the recursive call e is added to tt s . Given that the tree edge incident 
to u is the last one to be returned by choose(C, u), there is no path that does not use this edge, 
therefore it is not necessary to call line [8] for this edge. □ 

A. 2 Lemma [6] 

Proof: We proceed in order as follows. 

1. We only output a solution in a leaf and we only do recursive calls that lead us to a solution. 
Moreover every node partitions the set of solutions in the ones that use an edge and the 
ones that do not use it. This guarantees that the leaves in the left subtree of the node 
corresponding to the recursive call and the leaves in the right subtree do not intersect. This 
implies that different leaves correspond to different paths from s to t, and that for each path 
there is a corresponding leaf. 

2. Each left branch corresponds to the inclusion of an edge in the path ir. 

3. Since we are in a biconnected component, there is always a left branch. There can be no 
unary node as a right branch: indeed for any edge of B u t there exists always a path from s 
to t passing through that edge. Since the tree edge is always the last one to be chosen, unary 
nodes cannot correspond to back edges and binary nodes are always back edges. 

4. From point [1] and from the fact that the recursion tree is a binary tree. (In any binary tree, 
the number of binary nodes is equal to the number of leaves minus 1.) 



A. 3 Lemma [TT] 

Proof: Consider the following partition Vx = {r} U V<l U V3 where: r is the root; V% is the set of 
vertices with degree 2 and; V3, the vertices with degree > 3. Since Hx is compacted DFS tree of 
a biconnected graph, we have that V2 is a subset of the leaves and V3 contains the set of internal 
nodes (except r). There are no vertices with degree 1 and d(r) > 2. Let x = ^2 v& v 3 ^( v ) anc ^ 
y = J2veV2 We can write the density as a function of x and y: 



□ 



\Ex\ 



x + y + d{r) 



\V X 



2{\V z \ + \V 2 \ + l) 
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Note that | V3 1 < | as the vertices in V3 have at least degree 3, | Va | = | as vertices in V2 have 
degree exactly 2. Since d(r) > 2, we derive the following bound: 

Iffy I > x + y + 2 
\V X \ ~ |x + y + 2 

Consider any graph with \Vx\ > 3 and its DFS tree rooted at r. Note that: (i) there are no 
tree edges between any two leaves, (ii) every node in V2 is a leaf and (iii) no leaf is a child of r. 
Therefore, every tree edge incident in a vertex of V2 is also incident in a vertex of V3. Since exactly 
half the incident edges to V2 are tree edges (the other half are back edges) we get that y <2x. 

With \Vx\ > 3 there exists at least one internal node in the DFS tree and therefore x > 3. 

. . . x + y + 2 
minimize 

p + y + 2 

subject to < y < 2x, 
x>3. 

Since for any x the function is minimized by the maximum y s.t. y < 2x and for any y by the 
minimum x, we get: 

\E X \ > 9x + 6 > 11 
\V X \ ~ 8x + 6 ~ 10 

□ 

A. 4 Lemma [181 

Proof: We use standard data structures (i.e. linked lists) for the representation of certificate C. 
Persistent versions of these data structures exist that maintain a stack of modifications applied to 
them and that can restore its contents to their previous states. Given the modifications in /, these 
data structures take 0(|/|) time to restore the previous version of C. 

Let us consider the case of performing lef t_update(C, e). We cut at most 0(|Vx|) edges from 
C. Note that, although we conceptually remove whole branches of the DFS tree, we only remove 
edges that attach those branches to the DFS tree. The other vertices and edges are left in the 
certificate but, as they no longer remain attached to B u j, they will never be reached or explored. 
In the case of right_update(C, e), we have a similar situation, with at most 0(1^x1) edges being 
modified along the spine of the recursion tree. □ 
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(a) Bead string B Ut t at the root of the spine (b) Spine of the recursion tree 



Figure 4: Example bead string and spine of the recursion tree 

B Extended analysis of operations in a spine of the recursion tree 

In this appendix, we present all details and illustrate with figures the operations right_update(C, e) 
and lef t_update(C, e) that are performed along a spine of the recursion tree. In order to better 
detail the procedures in Lemma [TBI and Lemma [T7l we divide them in smaller parts. Fig. |4] shows 
(a) an example of a bead string B u ^ at the first node of the spine and (b) the nodes of the spine. 
This spine contains four binary nodes corresponding to the back edges in lb(u) and an unary node 
corresponding to the tree edge (u, v). Note that edges are taken in reverse DFS order as defined in 
operation choose(C, u). 

As a consequence of LemmaHl the impact of operations right_update(C, e) and lef t_update(C, e) 
in the certificate is restricted to the biconnected component of u. Thus we mainly focus on main- 
taining the compacted head Hx = (Vx, Ex) of the bead string B u>t . 

B.l Operation right_update(C, e) in a spine of the recursion tree 

Let us now prove and illustrate the following lemma: 

Lemma 20 (Lemma [16j restated) In a spine of the recursion tree, operations right_update(C, e) 

can be implemented in 0(\Vx\) total time. 

In the right branches along a spine, we remove all back edges in lb(u). This is done by starting 
from the last edge in lb(u), i.e. proceeding in reverse DFS order. In the example from Fig. HJ we 
remove the back edges (zi,u), (z2,u), (z^,u) and (z^,u). To update the certificate corresponding 
to B u t, we have to (i) update the lowpoints in each vertex of Hx (ii) prune vertices that are no 



14 



longer in B u j after removing a back edge. Note that, for a vertex w in the tree, there is no need 
to update 7(10). 

Let us consider the update of lowpoints in the DFS tree. For a back edge bi = (,2j, u), we traverse 
the vertices in the path from Z{ towards the root u. By definition of lowpoint, these are the only 
lowpoints that can change. Suppose that we remove back edge (24, u) in the example from Fig. H] 
only the lowpoints of the vertices in the path from 24 towards the root u change. Furthermore, 
consider a vertex w in the tree that is an ancestor of at least two endpoints Zi^Zj of back edges 
bi, bj. The lowpoint of w does not change when we remove b%. These observations lead us to the 
following lemma. 

Lemma 21 In a spine of the recursion tree, the update of lowpoints in the certificate by operation 
right_update(C, e) can be done in 0(\Vx\) total time. 

Proof: Take each back edge bi = (zi,u) in the order defined by choose(C, u). Remove bi from lb(u) 
and ab(zi). Starting from z%, consider each vertex w in the path from z% towards the root u. On 
vertex w, we update lowpoint(w) using the standard procedure: take the endpoint y of the first edge 
in ab(w) (the back edge that goes the nearest to the root of the tree) and choosing the minimum 
between 7(2/) and the lowpoint of each child of w. When the updated lowpoint (w) = 7(u), we stop 
examining the path from z^ to u since it implies that the lowpoint of the vertex can not be further 
reduced (i.e. w is both an ancestor to both Zi and Zj+i). 

The total cost of updating the lowpoints when removing all back edges bi is 0(|VxD : there are 
0(| Vx|) tree edges and we do not traverse the same tree edge twice since the process described 
stops at the first common ancestor of endpoints of back edges bi and By contradiction: if a 

tree edge (x,y) would be traversed twice when removing back edges bi and 6j+i, it would imply 
that both x and y are ancestors of Zi and Zi+\ (as edge (x, y) is both in the path Zi to u and the 
path Zi+i to u) but we stop at the first ancestor of both z% and ^+1. □ 

Let us now consider the removal from the certificate of vertices that are no longer in B u j as 
consequence of operation right_update(0, e) in a spine of the recursion tree. By removing a back 
edge bi = (zi,u), it is possible that a vertex w previously in Hx is no longer in the bead string B u j 
(e.g. w is no longer biconnected to u and thus there is no simple path u ~» w t). 

Lemma 22 In a spine of the recursion tree, the branches of the DFS that are no longer in B U) t 
due to operation right_update(C, e) can be removed from the certificate in 0(\Vx\) total time. 

Proof : To prune the branches of the DFS tree that are no longer in Hx , consider again each vertex 
w in the path from z% towards the root u and the vertex y, parent of w. It is easy to check if w is 
an articulation point by checking if the updated lowpoint (w) > 7(2/). If that is the case and w is 
not in the leftmost path of the DFS, it implies that w £ B u t, and therefore we cut the subtree of w 
and bookkeep it in / to restore later. Like in the update the lowpoints, we stop examining the path 
Zi towards u in a vertex w when lowpoint (w) = 7(1*) (the lowpoints and biconnected components 
in the path from w to u do not change). When cutting the subtree of w, note that there are no 
back edges connecting it to B Ut t (w is an articulation point) and therefore there are no updates to 
the lists lb and oh of the vertices in B u ^. Like in the case of updating the lowpoints, we do not 
traverse the same tree edge twice (we use the same halting criterion). □ 

With Lemma [21] and Lemma [22] we finalize the proof of Lemma [16] Fig. [5] shows the changes 
the bead string B u t from Fig. [4] goes through in the corresponding spine of the recursion tree. 
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(a) B Uit after removal of (zi,u) (b) B u , t after removal of (zi,u), (z2,u) 




(c) B Uyt after removal of (zi, u), (z2,u), (zs, u) (d) B„ it at the last node of the spine 



Figure 5: Example application of right_update(C, e) on a spine of the recursion tree 
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Figure 6: Block tree after removing vertex u 



B.2 Operation lef t_update(C, e) in a spine of the recursion tree 

To compute lef t .update (C, e) in the binary nodes of a spine, we use the fact that in every left 
branching from that spine, the graph is the same (in a spine we only remove edges incident to u and 
on a left branch from the spine we remove the node u) and therefore its block tree is also the same. 
In Fig. [H we show the resulting block tree of the graph from Fig. [4] after having removed vertex 
u. However, the certificates on these left branches are not the same, as they are rooted at different 
vertices. In the example we must compute the certificates C\ . . . C4 corresponding to bead strings 
B Zl)t ■ ■ ■ B Z4)t - We do not take into account the cost of the left branch on the last node of spine 
(corresponding to B v t) as the node is unary and we have shown in Lemma [T5l how to maintain the 
certificate in O(l) time. 

By using the reverse DFS order of the back edges, we are able to traverse each edge in Hx only 
an amortized constant number of times in the spine. 

Lemma 23 (Lemma [17] restated) The calls to operation lef t_update(C, e) in a spine of the re- 
cursion tree can be charged with a time cost of 0(\Ex\) to that spine. 

To achieve this time cost, for each back edge hi = (zi, u), we compute the certificate correspond- 
ing to B Zut based on the certificate of B Zi _ u t. Consider the compacted head Hx = (Vx,Ex) of 
the bead string B U)t . We use 0(\Ex\) time to compute the first certificate C% that corresponds to 
bead string B Zlt t- Fig. [7] shows bead string B Zl j from the example of Fig. |H 

Lemma 24 The certificate C\, corresponding to bead string B zi t, can be computed in 0{\Ex\) 
time. 

Proof: Let t' be the last vertex in the path u ~^ t s.t. t' E Vx- Since t' is an articulation point, the 
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subtree of the DFS tree rooted in t' is maintained in the case of removal vertex u. Therefore the 
only modifications of the DFS tree occur in head Hx of B U) t- 

To compute C\, we remove u and rebuild the certificate starting form z\ using the algorithm 
from Lemma [T3l restricted to Hx and using t' as target and j(t') as a baseline to 7 (instead of the 
depth). In particular we do the following. In order to set t 1 to be in the leftmost path, we perform 
a DFS traversal of graph Hx starting from z\ and stop when we reach vertex t'. We then compute 
the DFS tree, traversing the path z\ ~~> t' first. 

Update of 7. For each tree edge (v,w) in the t' z\ path, we set y(v) = ~y(w) — 1, using j(t') 
as a baseline. During the rest of the traversal, when visiting vertex v, let w be the parent of v in 
the DFS tree. We set 7(1?) = "f(w) + 1. This maintains the property that 7(f) > 7(10) for any w 
ancestor of v. 

Update of lowpoints and pruning of the tree. Bottom- up in the DFS-tree of Hx, we compute the 
lowpoints using the lowpoints of the children. Let z be the parent of v. If lowpoint{v) > lowpoint(z) 
and v is not in the leftmost path in the DFS, we cut the subtree of v as it does not belong to B zx ^. 

Computing lb and ab. During the traversal, when finding a back edge e = (v,w), if w is a 
descendant of v we append e to both lb(v) and ab(w); else we append e to both ab(v) and lb(w). 
This maintains the DFS order in the back edge lists. 

This procedure takes 0(|i?x|) time. □ 

To compute each certificate Cj, corresponding to bead string B Zi) t, we are able to avoid visiting 
most of the edges that belong B Zi _ t t- Since we take Zi in reverse DFS order, on the spine of the 
recursion we visit 0(|i?x|) edges plus a term that can be amortized. 

Lemma 25 For each back edge bi = (zi,u) with i > 1, let Ex[ be the edges in the first bead in 
common between B Zit t and B Zi _ x f The total cost of computing all certificates B Zi t in a spine of 
the recursion tree is: 0{Ex + Yli>i Ex'i)- 

Proof : Let us compute the certificate Cj : this is the certificate of the left branch of the iih. node of 
the spine where we augment the path with the back edge bi = (zi,u) of lb(u). 

For the general case of with % > 1 we also rebuild (part) of the certificate starting from z% 
using the procedure from Lemma [13] but we use information gathered in Cj_i to avoid exploring 
useless branches of the DFS tree. The key point is that, when we reach the first bead in common 
to both B Zi j and B Zi _ lt t, we only explore edges internal to this bead. If an edge e that leaves the 
bead leads to t, we can reuse a subtree of Cj_x- If e does not lead to t, then it has already been 
explored (and cut) in Cj_i and there is no need to explore it again since it is going to be discarded. 

In detail, we start computing a DFS from Z{ in B u ^ until we reach a node t' € B Zi _ x i. Note 
that the bead of t' has one entry point and one exit point in Cj_x- After reaching t! we proceed 
with the traversal using only edges already in Cj_i. When arriving at a vertex w that is not in the 
same bead of t', we stop the traversal. If w is in a bead towards t, we reuse the subtree of w and 
use 7(w) as a baseline of the numbering 7. Otherwise w is in a bead towards and we cut this 
branch of the certificate. When all edges in the bead of t' are traversed, we proceed with visit in 
the standard way. 

Given the order we take bi, each bead is not added more than once to a certificate Cj, therefore 
the total cost over the spine is 0(\Ex\)- Nevertheless, the internal edges Ex[ of the first bead in 
common between B Zit t and B Zi _ lt t are explored for each back edge bi. □ 

Although the edges in Ex\ are in a common bead between B Zut and B Zi _ 1; t, these edges must 
be visited. Since the entry point in the common bead can be different for Zi and £j-i, the DFS 
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tree of that bead can also be different. For an example, consider the case where Zi, . . . ,Zj are all in 
the same bead after the removal of u. The bead strings B Zi j . . . B Zj j are the same, but the roots 
Zi, . . . ,Zj of the certificate are different, so we have to compute the corresponding DFS of the first 
bead \j — i\ times. Note that this is not the case for the other beads in common: the entry point 
is always the same. 

Lemma 26 The cost 0(Ex + Si>i Bx'i) on a spine of the recursion tree can be amortized to 
O(Ex) in that spine. 

Proof: We can charge the cost 0(\Exi\) of exploring the edges in the first bead in common between 
B Zit t and B Zi _ lt t to another node in the recursion tree. Since this common bead is the head of at 
least one certificate in the recursion subtree of the left child of the ith node of the spine. Specifically, 
we charge the first and only node in the leftmost path of the ith child of the spine that has exactly 
the edges Ex'i as head of its bead string: (i) if \Ex[\ < 1 it corresponds to a unary node or a leaf 
in the recursion tree and therefore we can charge it with 0(1) cost; (ii) otherwise it corresponds 
to a first node of a spine and therefore we can also charge it with 0(\Exi\)- We use this charging 
scheme when i 7^ 1 and the cost is always charged in the leftmost recursion path of ith node of the 
spine, consequently we never charge a node in the recursion tree more than once. □ 

Lemmas 1251 and 1261 finalize the proof of Lemma [T7l Fig. [7] shows the certificates of bead strings 
B Zit t on the left branches of the spine from Fig. HJ 
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(d) Certificate of bead string B Zitt 
Figure 7: Certificates of the left branches of a spine 
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